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1. Introduction 

Lattices are basic algebra that provide a natural way to formalize and study the ordering 
of objects. They play an important role in various fields such as mathematics, computer 
science and engineering [1]. For example, they have applications in combinatorics, 
number theory, group theory, graph theory, set theory and logic circuits. 

Lattices have binary operations on elements of (partially) ordered sets called join 
(V) and meet (A). By these operations, we can investigate structure of ordered sets and 
manipulate elements of them algebraically. An application of Boolean lattice to logic 
circuits is one of successful case of such manipulation. 

Meanwhile, there have been numbers of researches on complex dynamical systems 
that are described by differential equations or difference equations. Most of complex 
dynamical systems exhibit either stable behavior, periodic behavior, or unpredictable 
behavior called chaos, depending on initial values and parameter of the systems. To 
classify the behavior of the systems, several entropy-based complexity measures have 
been proposed and employedjH [3]. 

In this article, we propose dynamical systems defined on algebra of lattices, which 
we call 'lattice equations'. We give exact general solutions of initial value problems 
for a class of lattice equations, and evaluate the complexity of the solutions. Since 
lattice equations are reduced to a certain class of max-plus equations or binary cellular 
automata (CA)|ll[5l|6] when we restrict dependent variables to the ordered set (M, <) or 
binary numbers, exact solutions of them apply to max-plus equations and CA. Though 
numerical or statistical studies about the behavior of solutions have been vastly done 
for max-plus equations and CA, direct evaluation of behavior of solution by means of 
explicit expressions is quite few. Therefore, our approach provides a novel viewpoint 
to analysis of complex dynamical systems, that is, analytical mathematics on lattice 
equations. 

We note our result is related to ultradiscretization procedure in the field of 
integrable systems that connects discrete max-plus equations, CA, and continuous 
differential or differences equations[71|8l|9l[T0]. The solutions are also directly connected 
one another through the procedure, though it is mainly found for the integrable systems. 
In this context, the lattice equations proposed with exact solutions will be useful in 
finding wider perspective between continuous and digital systems. 

The contents of this article are as follows. In Section |2l we first introduce the 
general definition of lattice and define a class where the complexity of their solutions 
are of polynomial order. Secondly we introduce the definition of elementary cellular 
automata (ECA) and show the relation between lattice equations and ECA. In Section |3l 
we show a list of lattice equations of the above class together with their exact solutions. 
In Section Hj we show some examples analyzing the asymptotic behavior of ECA using 



the exact solutions. In Section |5l we give concluding remarks. In Appendix A we give 



a list of the lattice equations and corresponding ECA which are equivalent to each other 



under the transformation of variables. In Appendix B we give some proofs about the 
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2. Lattice and Elementary Cellular Automaton 
2.1. Lattice 

A partially ordered set (poset) (L, <) is a lattice if the supremum (join) and the infimum 
(meet) of {a, b} always exist for any a, b E L[T]. Let us define V and A by 



a V 6 = sup{a, 6}, a Ab = mi{a,b}. (1) 
We have the following laws for any lattice, 

a\/ b = by a, a A b = b A a, (commutative laws) (2) 

a\/ {by c) = {ay b)\/ c, (3) 

a A {b A c) = {a A b) A c, (associative laws) (4) 

a V (a A 6) = a, a A (a V 6) = a, (absorption laws) (5) 

ay a = a, a A a = a. (idempotent laws) (6) 

Note that we can define V and A by 

a V 6 = max(a, 6), a A 6 = min(a, 6), (7) 
for the totally ordered set. 



The poset {L, <) is a 'distributive lattice' if it satisfies (one of) the following 
equivalent distributive laws for any a, b and c E L, 

ay {b Ac) = {ay b) A {aye), 

a A {by c) = {a Ab) y {a A c), (distributive laws) (8) 

{ay b) A {by c) A {cy a) = {a Ab) y {b A c) y {c A a). 

For example, the totally ordered set (M, <) for the real number field M with the usual 
order < is a distributive lattice since the above distributive laws hold. 

Let us consider a 'conjugate' element a for each a G L by the conditions, 

a & L, a = a, a >b if a < b. (9) 

Note that conjugate elements can not be defined for a general lattice and the definition 
is not unique even if they can be. For example, we can define the conjugate element for 
(M, <) by 

a = c — a (c: constant) (10) 

and also by 

{— ca (a < 0) 

-a/c {a > 0) 
where c > 0. For the conjugation, the following laws always hold, 

a V 6 = a A 6, a A 6 = a V 6. (12) 
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2.2. Evolution equations expressed by lattice operations 

Let us consider the following form of evolution equations using lattice operations, 

u]^' = f{u]_„u],u]^,), (13) 

where j denotes an integer space site number, n an integer time and u E L. The function 
/(a, 6, c) is constructed by V, A and ~. We call the equation in the form of f|T3|) 'lattice 
equation' for short. Let us consider the initial value problem of (1131) from the initial 
data assuming n = is an initial time without loss of generality. If we use the 
above equation recursively, we obtain the formal solution expressed by u° as 

where /*•"'•* is defined recursively by 

= /(/('=)(a_fc„i, . . . , afc_i), f^\a^k, . . . , a^), f^^\a_k+u • • • , dk+i))- (15) 

However this solution contains 3" terms of m° and it does not give any information on 
the solution since it is formal. In this article, we show the list of evolution equations 
with a general solution containing the terms of which the number is of polynomial 
order. If the solution can be expressed by the 0{n^) terms among {u^}, let us 
classify as the equation is of class Pm (polynomial class of order m). 
For example, let us consider the equation 

Using the basic laws of lattice shown above, we can obtain a reduced expression of 
solution, 

< = A A ... A A m5+„. (17) 

Equation f[T^ is of class Pi since the RHS of f[T7j) contains 2n + 1 (= 0{n)) terms which 
is equal to the infimum of the terms from to 

There are equivalent equations through the transformation of variables and 
coordinates. If /i(a, 6, c) = /2(c, 6, a), then the solution to n!^^^ = fi{u^_i, , u'^_^_i} 
is obtained from the solution to m""*"^ = /2(m"_i, u", m"+i) through the transformation 
j —7- —j. Similarly, the following evolution rules given by /j's are equivalent each other. 



fi{a,b,c) = f2{c,b,a), (reflection) (18) 

fi{a,b,c) = f2(a,b,c), (conjugation) (19) 

/i(a, *, *) = /2(*, b, *) = fs^*, *, c), (Galilean transformation) (20) 

/i(a, 6, *) = /2(*, 6, c), (Galilean transformation) (21) 

/i(a, *, c) = /2(a, b, *), (separation of even and odd sites) (22) 

where the symbol of arguments denotes fi does not depend on the corresponding 
argument. Moreover if /i and /2 satisfy the condition 

fi{a,b,c) = f2{a,b,c) = /2(a,6,c), (23) 
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they are equivalent through the transformation — ^27+1 ~^ '^2j+v Similarly, 
if /i and /2 satisfy the condition 

fi{a,b,c) = /2(a,6,c) = /2(a,6,c), (24) 

they are equivalent through the transformation m^" — )• m^" and -u^' 



2.3. Elementary cellular automaton 

Elementary cellular automaton (EGA) is a simple evolutional digital system and the 
evolution rule is given by 

uf' = Mu^^lu],u]^,), (25) 

where j denotes an integer space site number, n an integer time and u a binary state 
value (= or l)[5l[6]. The binary- valued function ff, determines the time evolution rule 
and gives the state value at the next time. Thus fb can be defined by a rule table in the 
following form. 



aba 


111 


110 


101 


100 


oil 


010 


001 


000 


fb{a, 6, c) 


bj 


be 


&5 


&4 


b3 


&2 


bi 


bo 



(26) 



where the upper row shows combinations of binary variables a, b, c and the lower the 
binary values of /(a, b, a). Since / is uniquely defined if all b^s are given, we can identify 
the evolution rule by the rule number (r)io = (&7&6 • • • ^1)2 after Wolfram [5l |6]. Let us 
call the EGA of the decimal rule number r 'EGAr'. There exist 256 rules from EGAO 
to EGA255. 

Let us consider the lattice (M, <) and assume that the conjugation operator is given 

by 

a = 1 - a. (27) 

If we define the rule fb{a,b,c) using V, A and ~, they are equivalent to the bitwise 
operations AND, OR and NOT respectively for the binary values as follows: 

a V 6 = max(a, b) = a OR b, a Ab = min(a, b) = a AND 6, 

a=l-a = NOTa, (28) 

where a,b & {0, 1}. Therefore every EGA rule can be expressed by the lattice operations. 
Moreover, if we consider the lattice (M, <), flT6l) becomes 



(29) 



and it gives EGA128 in the binary case. The solution 

,0 „,o 



u 



is 1 if 



u 







min(M°_„, 



5 '^j+n-li "^j+nJi 



(30) 



u 





j+n-l 



U 





j+n 



1 and otherwise in the binary case. 
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Considering the equivalent relations described in Subsection \'2.'2\ there exist 75 
independent rules of EGA. We list the rule numbers as follows: 

0, 1, 2, 3, 4, 6, 8, 9, 11, 12, 13, 14, 18, 19, 21, 22, 24, 25, 26, 

27, 28, 29, 30, 32, 33, 35, 36, 37, 38, 40, 41, 42, 43, 44, 45, 

46, 50, 54, 56, 57, 58, 60, 62, 72, 73, 74, 76, 78, 94, 104, 106, (31) 

108, 110, 122, 126, 128, 130, 132, 134, 136, 138, 140, 146, 

150, 152, 154, 156, 162, 164, 168, 170, 172, 184, 200, 232. 

3. Solutions to evolution equations expressed by lattice operators 

In this section, we list up the lattice equations satisfying the following conditions, 
(a) Equation form obeys ( IT3ll . 



(b) The terms a, b and c do not appear or appear once in /(a, b, c). 

(c) The solution is of polynomial class. 

In addition, we show some equations which satisfy (a) and (c) but not (b) as an 
exceptional case. All solutions are derived exactly by evaluating the initial value problem 
using lattice formulas shown in Subsection l2.1[ Two examples of proofs on such solutions 



are shown in Appendix Appendix B In the list blow, we show the form of /(a, b, c), the 



rule number of corresponding EGA and its solution. 

3.1. class Pq 

• /(a,6,c) = e (G L) (EGAO) 

• /(a, 6, c) = a (EGA240) 

J J-n- 

• f{a,b,c) = aAb (EGA12) 



u] = u'^_, A u]. 
/(a,6,c) = a A6 (EGA3) 

u] = A ^, = V (u°_2 A mJ) , 

f{a,b,c) = aAbAc (EGA2) 



/(a,6,c) = a A6Ac (EGA4) 



u'f = A M° A ul 



+1- 
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2m-l 



U 



2m 



(^i°_2 V V «°) A V «° V «° i) A («; v ^ v m° 2)- 



f{a,b,c) = (a V 6) A c (ECA42) 



f{a, b,c)^{aVc)Ab (ECA200) 

< = K_iVx.;+ja4 

f{a,b,c) = (a Vc) A6 (ECA76) 

= (^v^)a4 

/(a, 6, c) = (a V c) A 6 (ECA19) 



/(a, 6, c) = (6 V (a A c)) A (6 V (a A c)) (ECA36) 



u] 



V A ul,)) A V A «^), 

(«? V^?-! A u'^^ A K V A ^) 

A V V ^ V (xi;_2 A 

A V ^ V lij+i V A 4^)) 

AK_iV<:;vii;+iV^)) (n>2). 



5.^. class Pi 

We define new symbols about V and A as 

V Xk = Xfc,, V Xko+i V • • ■ V Xfcj , 

ko<k<ki 

A Xfc = Xfeo AXfco+i A • • • AXfci- 
ko<k<ki 

• f(a,b,c) = a A6 (ECA192) 

= A M° , 

= M°_„ A A . . . A mJ. 

• /(a, 6, c) = a A 6 A c (ECA128) 

li" = A 

~ ""j-n A itj-^+i A ... A u^_^_^_l A 
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f{a,b,c)^aAbAc (ECA8) 



u':={ A <+ J A A M° 

-l<k<n-2 ^^^) \0<fc<n 



• f{a, 6, c) = a A 6 A c (ECA32) 



= A mo_„+i a m;_„+2 a m°_„+3 a ... a a 



5.5'. class P2 
• /(a, 6, c) = (a V 6) A c (ECA138) 



0<l<n\\l-l<k<n-2 ^^'^J 



A(^V^V...V«o+„_2V«;+J 



A (m°+„_3 V m0+,_2 V M?+„_2) 
A (li°+„_2 V li.V-l) A U^n- 

f(a, 6, c) = (a V 6) A c (ECA162) 
aK -n+2 

Vi.0 

-n+3 -n+5 
A^V4V <n-3 V^g+n-l) 

A (ii°+„_2 V A 

/(a, 6, c) = (a V c) A 6 (ECA140) 

= (i/0_i V V ... V iiO+„_3 V ii0+„_2 V 1^;+ J 
A V ^ V ... V «o+„_3 V 



A(«o_iVmOVm;+2; 



A(mO_iVmO+i) a4 
/(a, b,c)^{aVc)Ab (ECA50) 



^ lo<;<m-lV ^ ^rn+zi+L V-m+;+l<A;<m-;-l ^^^^J 
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~ (''^j-2m+l ^ '^^-2m+2 ^ '"i-2m+4 A ... A 'Uj+2m-4 A 'W°+2m-2) 

V ('"j-2m+2 A 1*j-2m+4 A ... A 'U°^2m-4 A 1*j+2m-2 A Mj,+2m-l) 
('"'j-2m+2 A tij_2m+3 A 'iij-2m+5 A ... A W^+2m-5 A 1ij+2m-3) 

^ ('"'j-2m+3 ^ '^j-2m+5 A ... A U^^2m-5 ^ '^j+2m-3 ^ '"'j+2m-2) 

V A A lij+i) V A lij+i A ^) 
V^°_i A^) V(^A?i°+i), 

{o<y<. ( ^) ^ } 

{o<K^.-l (^^^^ ^ (-.+.-,l4,.-.-l "^^^) ) } 

{'^j-2m A Mj_2m+1 A W°_2m+3 A ... A 'U^+2m-3 A U^jj^2m-l) 
^ (""^-2™+! A 'Uj_2m+3 A ... A W°^2m-3 ^ ""^+2^1-1 A '"j'+2m) 

V (Mj_2m+1 A 'U^_2^_|_2 A 'Wj_2m+4 A ... A 'U°^2m-4 A 1lj+2m-2) 
('"'j-2m+2 A U^j-2m+i A ... A tij+2m-4 A U^+2m-2 A '"'j+2m-l) 

V (ii5_2 A ^ A ^) V A ^ A ii°+2) 

v(^ax.;)vKa^). 



V 
V 

V 



• /(a, 6, c) = (a V 6) A c (ECA21) 



A<1 A <,,_iVliO V V <+2fe+i+3 

.0<Km-2V ^ V0<fc<m-i-2 



A 



{o<l<m-2 ( (o<fe<m-i-2 ^^2^+^+2) V «°+2m-/-l) } 

v(u^,A.^,a(^^^A_^«;,,.2 
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^ {o<,X.-(C,.<^.-.- "?«".-')} 
• /(a, 6, c) = (o V 6) A c (ECA69) 



„2m _ /.,0 



(u»Au»^l) 

v{»»_,a<:a(^^a^<„)} 

/(a, b, c) = (a V 6) A (6 V c) A (c V a) (ECA232) 

= V A V 
A («;_2 V V u%,) A V uj+i V m;+2) 
A («5_3 V «;_2 V «5 V «5+2) A («;_2 v «° v v «?+3) 

A V V tz°_„_,3 V ... V 

A V ^i;_„+3 V ... V u%^_, V J 
A V V ... V m5+„_2 V 
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• /(a, 6, c) = 6 A (a V c) A (a V c) (ECA132) 



A(mO_,Vm;+i) 
A(^V^V^Vm°+2) 

A V V ... V V < J 

A(^z;_2V^V^V^) 

A V n0_„+2 V m0_„+3 V ... V 



4. Asymptotic behavior of ECA 

For the solutions of polynomial class, we can easily grasp its asymptotic behavior for 
n ^ in the binary case. We show it using two examples in Section [31 

The evolution rule in the case of /(a, 6, c) = a A 6 A c is equivalent to ECA8 in the 
binary case. We show again the solution, 



u] = A u° A ... A M°+„_2 A u° A A ... A 
In the binary case, u Au = since m G {0, 1}. Thus the above solution is simplified as 

1 = and M° = = 1) 

(otherwise) 



u] = A M° A ul 



'i+i 



u] = {n> 2). 

Figure [1] shows a typical evolution pattern of ECA8. 

Another example is /(a, 6, c) = (a V c) A 6 (ECA140). The solution is 



A V ^ V ... V m0+„_3 V J 



(32) 



A (mO_, V V u%,) 
A(^Vm;+i) Am°. 

In the binary case, it is clear that = (n > 1) if = and that = 1 (ra > 1) 
if u'^^i = and m° = 1. Moreover, -u" = (n ^ 0) if there exist r > such that 
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n 
i 




Figure 1. Evolution pattern of ECA8. The symbols □ and ■ denote and 1 
respectively. 

u^^i '^%r — 1 u^+r+i — 0- Therefore, we can see the asymptotic 

behavior of for large enough n as 

fx (m • 1 = and m° = 1) 

< = i / , ' in » 0) 

I (otherwise) 

Note that m" = 1 (n > 1) if u° = 1 as a special case. Figure [2] shows a typical evolution 
pattern of ECA140. 



j 




Figure 2. Evolution pattern of ECA140. 

Let US consider the asymptotic behavior of solution from real valued initial data 
for the same /(a, b,c) = (a V c) A 6 of the lattice (M, <). The lattice operations a V 6, 
a Ab and a can be replaced by max(a, b), min(a, b) and 1 — a respectively for the lattice 
(M, <). Thus fl32|) can be rewritten as 

u] = min( max(l - 1 - m°, . . . , 1 - m°+„_3, 1 - m°+„_2, wj+n), 
. . . , 

max(l — Uj_i, 1 — Mj, . . . , 1 — 'Uj_|_j._2, 1 — 'Uj_(_^_]^, Mj_)_^_|_]^), 

max(l - 1 - . . . , 1 - M°+^_2, U^j+r)^ (33) 
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max(l - Uj_^, 1 - Uj,Uj_^2), 
max(l - 

If there exist r > 1 such that u^^^. < 1/2, the max terms above the underhned terms in 
(!33|) are all greater than max(l — 1 — . . . , 1 — ^'j+r) since they include 

the same 1 — 1 ~ '^j'^ • • 1 ~ Mj'+r._2; and 1 — m^^,^. which is greater than 'U°_^_j,. They 
can be neglected since RHS of fl33|) is a minimum of all max terms. Thus we obtain 

u] = min(max(l - 1 - m°, . . . , 1 - m°+^_2, 1 - 

max(l - 1 - . . . , 1 - M°+^_2, 

• • • ) 

max(l - 1 - M°+2)5 
max(l - 

and m""*"^ = m". Figure [3] shows an example of this type of solution. Otherwise u1 > 1/2 
for any k > j — 1. In this case, 

max(l - 1 - . . . , 1 - u°+^_2, = «?+r 

and we have 

u] = min(u°, . . . , u°+„), 
from fl33l) . Then all u become the same value min u1 for n ^ under the periodic 

^ ^ l<k<K " 

boundary condition for the space cites with a finite period K. In any case, the 
asymptotic behavior of solution is static (m^^^ = m") for n ^ 0. 




Figure 3. Evolution pattern of real valued solution. 

The behaviors of asymptotic solution shown in Figs. [2] and |3] are similar each other 
in a sense that they become static. Though the binary solution is a special case of real 
valued solution, it implies that EGA can become a good test of the corresponding lattice 
equation. 
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5. Concluding remarks 

We proposed the equations constructed from lattice operations V, A and ~ of which 
solutions are expressed by the initial data of polynomial order. Since EGA are embedded 
in these equations as a special case, we can grasp the asymptotic behavior of EGA easily 
using the solutions. 

We proposed 24 equations and their solutions. Rule numbers of corresponding EGA 
are as follows: 

0, 1, 2, 3, 4, 8, 12, 19, 21, 32, 36, 42, 50, 69, 76, 128, 132, 138, 140, 162, 192, 
200, 232, 240. 

Therefore about one third of EGA listed in ( 13T]) are solved. EGA are roughly classified 
into 4 classes (Glass 1 to 4) by Wolfram. Almost EGA expressed by the lattice equation 
shown in this article are in Glass 1 or 2, and some are in Glass 3. This classification 
is done by the geometrical complexity of solution pattern. On the other hand, we 
classified the lattice equations from a viewpoint of algebraic complexity, that is, the 
order of initial data included in the solution. It is interesting that some rule numbers 
in the same Wolfram class are not in the same polynomial class. 

The lattice operations in the lattice (M, <) can be embedded in the max-plus algebra 
where max (min) is the addition, + is the production, and — is the subtraction [1]. The 
max-plus equation 

m""^^ = Uj + max(M"_i, 1 - m") - max(M", 1 - 

can not be expressed only by the lattice operations and is equivalent to EGA 184 in the 
binary case[9]. However, this equation is transformed into the lattice equation of Pi, 

/7+i = max(/;^i,/]^i), 

through the transformation 

There exist some max-plus equations which reduce to the lattice equation of polynomial 
class like the above example. It is a future problem to clarify the relation between the 
max-plus and the lattice expressions from the viewpoint of solvability of initial data 
problem. 

We analyzed the lattice equations constructed from a simple combination of lattice 
operations. Though they are simple, we can not confirm the polynomial class of some 
equations, for example, those defined by the following /(a, b, c), 

(a V c) A 6 (EGA35) , (a V 6) A c (EGA168) , 
(aV6)Ac (EGA84), (aV6)Ac (EGA81). 

They may not be in a polynomial class. It is another future problem to evaluate the 
complexity of general solution to the above equations. In addition, there exist many 
equations defined by the more complicated combination of lattice operations. These 
propose a vast target to be analyzed. 
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Appendix A. Equivalent lattice equations and ECA rule numbers 

In this appendix, we show in Table lAll the equivalent lattice equations of polynomial 
class in the form, 

and its corresponding rule number of ECA. The first column of each row shows equivalent 
/(a,6, c)'s and the first one is a representative. The second column shows the class of 
its general solution. The rule number of corresponding ECA is noted together with 
f{a,b,c). Note that classes of some lattice equations in the simple form described in 
Section [3] is not yet confirmed. We show them in Table IA2[ 



e e L (0, 255) 


Po 


a (240), b (204), c (170) , a (15), b (51), c (85) 


Po 


aAb (192), bAc (136), a Ac (160), 


Pi 


a V 6 (152), 6 V c (238), a V c (250) 




aAb (12), bAc (34), a V 6 (207), 6 V c, (187) 


Po 


aAb (48), bAc (68), a V 6 (243), 6 V c (221), 




a Ac (10), a Vc (175), a Ac (80), a Vc (245) 




aAb (3), bAc (17), a V 6 (63), 6 V c (119), 


Po 


a Ac (5), a V c (95) 




aAb Ac (128), aVfoVc (254) 


Pi 


a Ab Ac (8), a Ab Ac (64), a V 6 V c (239), a V 6 V c (253) 


Pi 


a A 6 A c (32), a V 6 V c (251) 


Pi 


aA6Ac(2),aA6Ac (16), a V 6 V c (191), a V V c (247) 


Po 


aAb Ac (4), a V6 Vc (223) 


Po 


a A 6 Ac (1), a V 6 V c (127) 


Po 


(a V c) A 6 (200), (a A c) V 6 (236) 


Po 


(a V c) A 6 (140), (a V c) A 6 (196), (a A c) V 6 (206), (a A c) V 6 (220) 


P2 


(a V c) A 6 (76), (a A c) V 6 (205) 


Po 


(a V c) A 6 (50), (a A c) V 6 (179) 


P2 


(a V c) A 6 (19), (a A c) V 6 (55) 


Pi 


(a V 6) A c (138), a A (6 V c) (208), (a A 6) V c (174), a V (6 A c) (244) 


P2 


(a V 6) A c (162), a A (6 V c) (176), (a A 6) V c (186), a V (6 A c) (242) 


P2 


(a V6) Ac (42), a A (6 Vc) (112), (a A6) Vc (171), a V (6 Ac) (241) 


Po 


(a V 6) A c (69), a A (6 V c) (13), (a A 6) V c (93), aV {bAc) (79) 


P2 


(a V 5) A c (21), a A (6 V c) (7), (a A 6) V c (87), a V (6 A c) (31) 


P2 


(a V 6) A (fe V c) A (c V a) (232), (a A 6) V (6 A c) V (c A a) (23) 


P2 


6 A (a V c) A (a V c) (132), 6 V (a A c) V (a A c) (123) 


P2 


(6 V (a A c)) A (6 V (a A c)) (36), (6 A (a V c)) V (6 A (a V c)) (219) 


Po 



Table Al. Equivalent lattice equations and their class. 
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(a V c) A 6 (35), (a V c) A 6 (49), (a A c) V 6 (59), (a A c) V 6 (115) 
(a V 6) A c (84), a A (6 V c) (14), (a A 6) V c (213), a V (6 A c) (143) 
(a V 6) A c (168), a A (6 V c) (224), (a A 6) V c (234), a V (6 A c) (248) 
(a V 6) A c (81), a A (6 V c) (11), (a A 6) V c (117), a V (6 A c) (47) 



Table A2. Equivalent lattice equations of unconfirmed class. 



Appendix B. Proof of solutions 

The solutions described in the section |3] are all derived only by using the formulas in 
the subsection 12.11 In this appendix, we show two examples of proofs of solutions. The 
solutions to the other equations can also be proved by the similar procedure. 
The first example is about the equation 

u]^' = ufr,Au]Au]^„ (B.l) 

which corresponds to ECA8. The solution is 



< = A A ... A A u', A n^^i A • • • A w^V- (B-2) 



It is of class Pi. For n = 1, it is easy to check wj = u^_^ A m° A u^+i from (IB.ip and 
(IB.2I) . If we substitute it into the RHS of (IB. II) . we obtain 



A w'; A w;Vi 



(m0_2 a ... a m0_^„_3 a a ... a mO+„_i) 
A (^ A . . . A A u^, A ... A 

A (^ A . . . A A ^ij+i A ... A uJ+^+i 

(m;_2 V ... V ?x;+„_3 V . . . V ^zO+„_i) 



A («°_i A ... A a a ... a m°+„+i) 

= ((m°_2 V . . . Vm;+„„3 V . . . V mO^.„„i) 

A^) A (^. . . AmO^„_iAm; a . . .Au%^^,) 
= A ^0 ... A a n° a ... a 

Thus the induction on the solution holds. 

The second example is about the equation 

u]+' = i^,Vu]^,)Aul (B.3) 

which corresponds to ECA140. The solution is 



(mO_, V V ... V m0+„_3 V V 
A(^V^V... VnO+„_3V^i^.+„_i, 







(B.4) 
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It is more complicated than that of the first example and is of class P2. For n = 1, it is 



easy to check Wj = {u^^^ V u^^^) A m° from (IB.3P and (IB.4p . We can rewrite w'^ 



and Wj+i as follows: 



,0 



. V m0+„_3 V a ... a V a 

[u-^, A ... A m;+„_3 a V ... V a ^) V^, 



A B 



i = V ... V V n°+J A ... A {u^^_, V nj+i) A 

^ V ( V . . . V V J A ... A V u]^^) A ) ) 



A 



D 



B C 




(mO V . . . V V A ... A V u^^^,) A nj+i 



A V ... V V < J A ... A V A mJ+i . 



c 

If we substitute (lB.4p into the RHS of (IB.Sp . we obtain 



(^I^ V w^+i) A = {{A V fi) V A C)) A ((5 V C) A D) 

= (A V V ^) A (A V 5 V C) A (5 V C) A D 
= (A V 5 V A (5 V A D 



Moreover, we have 

AyE= A ... A m;+„_3 A ^^,) V ... V (m°_2 A ^) 
^ V ' 

A 

V(^V...VMO+„_iVn°+„+i) 

V ' 

= ((m;_2 a ... a m5+„_3 a V 

V... V ((n°_2A^) V^) Vn;+„+i 

= 

Therefore, 

(^Vt^;+i) A^; 
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{By E) A{BV C) AD 



,0 

'j+n+l 



,0 

'j+n-2 



,n+l 



is derived and the induction on the solution holds. 
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